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Abstract 

This is a review of the cosmological bounds on B — L violating interactions, and the loopholes 
in the argument that gives these constraints. If one assumes that the baryon asymmetry we 
observe today was present above the electroweak phase transition in equilibrium with the non- 
perturbative B + L violating processes, then interactions that violate all three of {B/3 — Li} 
cannot simultaneously be in equilibrium. Otherwise the baryon asymmetry would be washed 
out. Therefore violation of at least one of the B/3 — Li must be small. This argument can be 
evaded by not having the observed baryon asymmetry present in the thermal bath (for instance, 
make it at the electroweak phase transition), by using a non-standard cosmological model, or 
possibly by some mass effects in models where the difference between lepton flavour asymmetries 
is conserved. 


1 introduction 

It is observed that at least up to the scale of our local galaxy cluster, the Universe we see is made 
of matter rather than anti-matter. It is difficult to build a cosmological model where the Universe 
contains equal numbers of baryons and anti-baryons separated on scales larger than galaxy clusters 
(There is a problem with causality [Q), so one usually assumes that the Universe as a whole contains 
a net excess of baryons over anti-baryons, and that at some time in the early history of the Universe, 
this asymmetry was generated. This is the Baryon Asymmetry of the Universe, or BAU (for reviews, 
see, e.g. [||, |, 0.) 

There are three ingredients required to create a baryon asymmetry [0 : baryon number violation, 
C and CP violation, and some out of equilibrium dynamics. The Standard Model (SM) does contain 
all these ingredients, but not in the right quantities to generate the BAU. 

The Standard Model violates C maximally (because it is chiral), and contains small amounts of 
CP violation in the CKM matrix (probably not enough to generate the BAU though). It has no 
perturbative baryon number violation. The lowest dimensional baryon number violating operator 
allowed by standard model gauge symmetries and involving SM fields is qqql [0, which is of dimension 
six. However the Standard Model does contain non-perturbative baryon number violation |0: the 
current associated with the quantum number B + L (baryon -|- lepton number) is anomalous, so 
certain gauge field configurations can induce the B + L violating vertex (qqqi)^ (see, for instance, [0 
for a review and references). This operator allows processes such as 

UgdifdfpCgSbS'ptgbbhfp > (f) 

where g,r,b are SU(3) colour labels, not summed, and the quarks are SU(2) doublet members. At 
zero temperature the rate for such B + L non-conserving processes is exponentially suppressed by 
the action of the weak gauge field configuration, so one does not expect to see low energy B + L 


violation. At finite temperature the rate is much less suppressed; the rate above the electroweak 
phase transition is 

Tb^l = ( 2 ) 

where k is a constant of order 1, and n is 4 or 5 [^. This is clearly in thermal equilibrium in the early 
Universe {T > H = the expansion rate of the Universe; if T > hf the timescale for the interaction 
is less than the age of the Universe.). After the phase transition, the rate is Boltzmann suppressed 
(r ~ exp{—mw/(Q!T)}), so rapidly drops out of equilbrium. This is potentially very interesting for 
baryogenesis (as was noted by H]): a source of baryon number violation that is efficient in the early 
Universe when the asymmetry needs to be generated, and then turns off at low energy in perturbative 
processes where we see no baryon number violation. 

Unfortunately the Standard Model lacks the third ingredient required to generate the baryon 
asymmetry. There is not a big enough departure from equilibrium at the electroweak phase transition; 
the B + L violating processes turn off gradually and no asymmetry in generated |^. The observation 
of the baryon asymmetry can therefore be interpreted as evidence for new physics, since it does not 
appear possible to generate an asymmetry within the Standard Model. 

The new physics can have two effects on the generation of the baryon asymmetry. It may simply 
add enough CP violation and make the phase transition sufficiently strongly hrst order that the 
baryon asymmetry is generated at the phase transition using the SM B + L violation. Alternatively, 
the beyond-the-Standard Model physics may involve new sources of B and/or L violation which 
generate the observed asymmetry before the electroweak phase transition. In the hrst case there are 
no cosmological bounds on B — L violation. So for most of this review, I will assume the second 
possibility: suppose that the baryon asymmetry we observe today was generated early in the history 
of the Universe and is present as an asymmetry in the number density of quarks and antiquarks in the 
thermal bath just above the electroweak phase transition. The anomalous B + L violating processes 
are also operating in this plasma, so if the asymmetry generated previously wishes to survive in 
their presence, it must carry at least one of S/3 — L^^ S/3 — L^ or S/3 — L^-. These are the 
three global quantum numbers exactly conserved in the Standard Model. If there are simultaneously 
interactions in equilibrium that violate all three of these quantum numbers, all asymmetries will be 
washed out. Since we are assuming that the baryon asymmetry we observe today was generated 
before the electroweak phase transition, this is undesirable. One can avoid this outcome by requiring 
that processes violating one of the {B/?> — Lj} be out of equilibrium: 


Ts/a-Li <-ff for one z 


( 3 ) 


where H is the expansion rate of the Universe, and i is one of e, /r or r. This condition should hold 
for all temperatures during which the baryon asymmetry is at risk: 100 GeV ^ T ^ 100 TeV. The 
reasons for this range of temperatures will be discussed in section 3. 

For a yukawa-type B — L violating interaction (eg R-parity violating trilinears in Super- 

symmetry), equation (§) becomes 


lO-^A^T < 


25T2 


( 4 ) 


rripi 


which gives a bound of order A < 10“^. I use H = lOT'^/rUpi in the SM, and H = 2bT‘^/mpi for 
SUSY particle content. One can estimate the rate associated with a majorana neutrino mass to be 
of order 




10 


< H 


( 5 ) 


which gives my < 30 keV in one generation. The numerical factors in these bounds will be discussed 
more carefully in section 2. Bounds can also be set on non-renormalisable R/3—Lj violating operators 
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of dimension D = 4 + n > 4; if all the coupling constants are absorbed into the mass scale M, the 
rate can be estimated as 


10 


- 3 - 


M2" 


< 


rripi 


( 6 ) 


For D = n + A = Q, this gives M > IO'^Tq^y GeV. It is clear that non-renormalisable interaction rates 
increase faster with temperature than the expansion rate, so are more likely to be in equilibrium at 
higher temperatures. The maximum temperature at which the bounds can be applied, which gives 
the best bounds (T ~ 100 TeV), will be discussed in section 3. 

These constraints are comparatively stringent, so it us useful to know exactly where they apply 
and how they can be avoided. These bounds apply to only one of the {B/3 — Lj}, so lepton number 
violation in one family and baryon number violation must be small. In section 2 I will review 
estimates of interaction rates for various coupling constants, and the bounds one can derive from 
these estimates. The most obvious loophole is to generate the asymmetry at the electroweak phase 
transition, or more precisely, to not have the baryon asymmetry present in the plasma just above 
the EPT. In section 3 I will show that the BAU only needs to be absent from the thermal soup 
for temperatures between ~ 100 TeV and the EPT, and discuss various ways of accomplishing this. 
Section 4 contains various other mechanisms for avoiding the bounds such as low Treheat inflationary 
models. There is a summary of the constraints and loopholes in section 5. The remaining part of 
the introduction is a short overview of previous work on this topic. 

It has been known for a long time that if the baryon asymmetry is generated by the out-of- 
equilibrium decay of heavy GUT particles, then baryon number violating processes mediated by 
these particles must be out of equilibrium after the asymmetry is generated |^. Otherwise they 


will wipe out the asymmetry generated in the decay. For the baryon asymmetry to survive in the 
presence of the B + L violating electroweak effects, it must be an asymmetry in i? — L, which 
makes it difficult (but not impossible, see section 4.1) to preserve the baryon asymmetry in B — L 
conserving GUTs. An attractive alternative, introduced by Fukugita and Yanagidajl^, is to generate 
a lepton asymmetry in the decay of heavy singlet neutrinos, and use the electroweak B + L violation 
to transform this into a baryon asymmetry. Lepton number violating interactions mediated by the 
heavy singlet neutrinos must be out of equilibrium after the lepton asymmetry is generated, to avoid 
washing out the asymmetry||13||. This bound on lepton number violating processes is not limited 
to leptogenesis scenarios; Barr and Nelson [|^, and [^, observed that the baryon asymmetry, the 
B + L violating processes, and interactions that violate the three {B/3 — Li} cannot simultaneously 
coexist in equilibrium, and used this to set bounds on majorana neutrino masses in one generation. 
This argument was then extended to R parity violating interactions in the Supersymmetric version 
of the Standard Model [^, |T^ and to generic B — L nonconserving interactions [l^. The constraints 
and loopholes were then discussed by many people |^, |^, ^ ^ 


2 rate estimates and constraints 

The argument that gives cosmological bounds of B — L violation goes as follows: assume that the 
observed baryon asymmetry was generated before the EPT. To survive in the thermal soup, it must 
carry a quantum number that is “effectively conserved”— i.e. conserved by the interactions that are 
in chemical equilibrium. The SM interactions are all in chemical equilibrium just above the EPT, 
and the three global quantum numbers that they conserve are the {B/3 — Li}. To protect the BAU, 
any additional interactions present must also “effectively conserve” at least one of the {B/3 — Li}, 
or equivalently, interactions violating at least one of the {B/3 — Li} must be out of equilibrium. 

In the hrst part of this section, I will estimate interaction rates associated with coupling constants 
of different mass dimensions. It follows from these estimates that all the SM interactions come into 
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equilibrium before the EPT. In the second part, I will discuss the equations of chemical equilibrium 
and charge conservation for the SM interactions. These are interesting because there are small 
“second order” terms proportional to lepton masses that can be used to protect an asymmetry. 
Finally in the third subsection I will use the rate estimates from section 2.1 to set bonnds on neutrino 
masses, i?-parity violating coupling constants, etc. 

As one can see from the estimated bonnds in eqnations (^), (|]) and (|]), the best constraints 
are on renormalisable conplings. With the SM particle content, there are no B — L non-conserving 
operators of dimension less than 5. However, in the snpersymmetric extension of the SM, when im¬ 
parity is not imposed, baryon and lepton nnmber violating masses and trilinear terms are possible. 
In this proceedings, I will discuss constraints on B — L violation with SM particle content at the 
electroweak scale, and also with SUSY particle content. 

The MSSM superpotential is 

W = fiH.H, + + hl^H.L.E;; . (7) 

The Lagrangian also contains kinetic terms, gange interactions, D-terms and soft SUSY breaking 
terms of the form 


soft masses + BhHiH 2 + Aim^QpUl + A^fHiQpD^ + A^iH^UE^. 


( 8 ) 


I am abusively using capital letters for both superhelds (as in eqn and scalar component helds (as 
in eqn |^). Quark generation indices are p,q,r,s... and lepton indices are Whether indices 

are np or down makes no difference. 

The above SUSY Lagrangian is with the symmetry /^-parity imposed. Other renormalisable 
interactions can be constrncted with the MSSM particle content, bnt these new interactions violate 
either B or L, and are often removed by imposing a symmetry like i?-parity |^. See for other 
possibilties. These are examples of the interactions I am interested in setting bonnds on. There are 
possible new snperpotential terms 


Wfi = AH2L, + X^^L^L^El + 

and new B or L violating soft terms 

soft masses mixing Ll and Hi 


( 9 ) 


( 10 ) 


I assnme that SUSY masses are of order 100 GeV ~ Tf.. The bonnds on trilinears are not 
substantially weakened by larger soft masses—if one sets the constraints at T ~ rnsusY ~ 1 TeV, 
they are weakened by a factor of abont 3. The effects of larger SUSY masses were discussed in 

Bii. 


2.1 rate estimates 

In this section I will ronghly estimate interaction rates for a selection of operators. See for instance 


171 CT for more thorongh discussions that disagree in details bnt all get approximately the same 
constraints. For a detailed calcnlation of the interaction rate associated with a specific conpling 


constant see, e.g. 


The rate for a boson to scatter off a fermion via the yukawa interaction |/0'^i'02 can be estimated 

for an introdnction) to be m 


using kinetic theory (see e.g. 


r(0V'i ^ -< crvn^^ 


>- 


4 X 


( 11 ) 
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where a is the T = 0 cross section, is the equilibrium thermal distribution of fermions 'ifj i, v is 
the relative velocity of the boson and the fermion (~ c = 1 ), and the brackets represent averaging 
the cross section over thermal distributions in momentum space. 

For a gauge boson, the matrix element, and therefore the numerical coefficient, will be different, 
but if one substitutes the gauge coupling g for the yukawa h in equation (pH]), it is clear that gauge 
interactions are in equilibrium up to temperatures of order ^ mcuT- 

One can similarly estimate the rate for (p'ljji —>■ 'ijji'j where 7 is an arbitrary gauge boson, to be of 
order 

r(0'0i —>• ~ 5 X 10~^y^aT . (12) 


The rate for a scalar (Higgs or spartner) to decay to two fermions via the interaction y(t)'>p'ip can 
be estimated as [O 

-^ 9 9 

oV mi 

1.4 X 10 “^^—^ . 

T 


r 


(13) 


rricj, is the decaying particles mass in the hnite temperature effective potential. Equation (|T3|) neglects 
hnal state phases space suppresion due to fermion masses, which can be signihcant if the mass of the 
scalar is of order gT. (The hnal state fermions can also have thermal masses of order gT.) 

To estimate a rate associated with a yukawa, one therefore has a choice between scattering and 
decay. The scattering rate has an extra power of a, but the decay may be suppressed by the small 
amount of hnal state phase space, depending on the (thermal) masses of the participating particles. 
I will use the decay rate estimate to set bounds on B — L violating couplings, because the decaying 
particle is a SUSY spartner with a zero temperature mass of order 100 GeV. 

In the SM case, y ~ \/2mf/v where v = 246 GeV is the Higgs vev, and is a fermion mass. 
(This is not entirely accurate for the quarks, because the hnite temperature mass eigenstates will 
not be the same as the zero temperature ones, but this should be a reasonable estimate). Taking 


H ~ lOT^/rripi, one hnds that the electron yukawa is in thermal equilibrium below 




TeV [0, the ur quark will be in chemical equilibrium below ~ 100 TeV, and so on. 


10 - 100 
(See 


24 


for a careful discussion of when the electron yukawa comes into equilibrium, taking into account the 
thermal masses in the decay rate.). This means that all SM interactions are in equilibrium just above 
the EPT. 

To estimate a rate associated with a majorana neutrino mass, I need to consider how it is gen¬ 
erated. If the zero temperature mass arises from the dimension 5 operator y‘^{iiHH)/M, where M 
is a heavy right-handed singlet neutrino mass, then at temperatures Tc < T M this operator 
generates two-fermion-two-Higgs scattering. The rate can be estimated |]T^ as in equation (pT|) 
to be 


T{iH ^ £if) ~ 4 X 10' 


M2 


(14) 


This is a conservative estimate of the lepton number violating rate, because the latter should also 
include the processes £i —»• HH, HH —>■ i£ and so on. 

Alternatively, if the majorana mass is present above the EPT as a mass, then it generates lepton 
number violating processes such as uW —>■ uW. Majorana mass insertions on the external neutrino 
legs contribute the lepton number violation. The mass corrections to the scattering process can be 
estimated from ([TTD as 

Y{uW 9W) ~ IQ-yT (15) 

and to a decay as 

r(IT ^ J,z.) ^ 3 X 10-Vt (16) 
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(To make this estimate, I took rri^{T) ~ /?> and rnl{T) ~ g‘^T‘^/S. For these masses, the decay 

W ^ uu is kinematically not allowed, so I ought not to use this estimate. It is nonetheless in the 
literature. ) 


2.2 chemical equilibrium equations 


There are two approaches in the literature to the equations of chemical equilibrium as applied to 
cosmological contraints on B — L violation. The more common one uses kinetic theory |^, , and 

is reliable insofar as it agrees with the more esotheric calculations that start from the free energy 
p6| . I will use something in between, that I think gives the right answers as long as one is not to 
near the EPT (See for an analysis that is valid for all temperatures.) 

If an interaction is in chemical equilibrium, for example the decay (j) —> then the sum of 

the chemical potentials of the participating particles is zero: 




( 17 ) 


Above Tc, all the gauge bosons have zero chemical potential |T^, |2^. From equation (^) it 
is clear that SM gauge interactions are in equiliibrum up to very high temperatures, so all members 
of a gauge multiplet have the same chemical potential. 

At T ~ Tc, all the SM processes are in equilibrium. If the chemical potential for a gauge multiplet 
is written as the name of the multiplet (so the chemical potential for is then the Yukawas 
imply 


—+ t + H = 0 
—u^ + — H = 0 

-d)^ + + H = 0 


(18) 


The rate estimates in section 2.1 say that the electron yukawa is in equilibrium below ~ 10 — 100 
TeV. So equations (|I^ apply for Tg < T < 10 — 100 TeV. Above this temperature, the cr have a 
thermal distribution (they have gauge interactions) but any asymmetry in the cr is decoupled from 
the rest of the plasma. In the absence of an interaction that can transfer the cr saymmetry to other 
particles, asymmetries must remain amoung the other particles (in particular, a baryon asymmetry) 
to cancel the hypercharge carried by the cr asymmetry. 

The non-perturbative B + L violating processes give 


9gL + E^* = o 

i 


(19) 


Note that these interactions eat three units of baryon number and three units of lepton number, so 
they violate B + L but not B — L. However, when they are in equilibrium with SM interactions in 
the early Universe, they imply ur = and not hr = —This is because the interactions 

in equilibrium try to minimise the free energy, and the minimum is not neccessarily where all the 
quantum numbers are zero. The “sphalerons” of equation (|^) do not take hr^r to zero because they 
only eat SU(2) doublets, and the SU(2) singlets (who also carry B or L) have different hypercharge. 
In a plasma with Y = 0 but non-zero ur-r, the free energy minimum is at ur = —(28/51)nL rather 
than at ur+r = 0. 

In the case of the MSSM, there are many new particles. The gauginos have majorana masses 
(sufficiently large to be in chemical equilibrium) so can carry no asymmetry and must have zero 
chemical potential. This means that the sfermion and fermion chemical potentials are equal and 
opposite to each other, because the gaugino interaction gcfy^ip imply 

0-FV’ + ^ = 0 + V’ = O (20) 
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The two Higgs have opposite chemical potentials due to the mass term so one gets the same 

chemical potentials as in the SM, subject to the same constraints due to yukawa couplings. 

If the flavour off-diagonal slepton masses are large enough for lepton flavour violation to be in 
equilibrium, then only B — L is conserved, and not the three H/3 — Lj separately. I will assume 
for the moment that the off-diagonal masses are “small enough” < 10“^), and come back 

to this constraint when I am discussing bounds on B — L violating operators in supersymmetric 
theories. 

This gives six equations for the 10 unknown chemical potentials. The four free parameters 
correspond to the four conserved charges: electric charge (or hypercharge) and the three {H/3 —Lj}. 
The electric charge must be zero (imposing H = 0 gives an equivalent constraint): 

^ Qemin - h) = 0 (21) 

particles 


where {n) n is the (anti-) particle number density, n — n = dF/dfi where F is the free energy and fi 
the particle chemical potential. One hnds 


‘^9 rp 2 

n — n = 


1 + 0{m‘^/T^) 

2 + Olm^/T^) 


fermions 

bosons 


( 22 ) 


for m T. g here is the number of degrees of freedom of the field—2 for a chiral fermion and a 
charged scalar. One can correctly get the first term of equation (^) using kinetic theory, by expanding 
the equilibrium particle distributions in fi and (see section 4). However, this expansion does not 
always give the right coefficient for the 0{fim^) terms. Note that I am working above the EPT, 
so the masses in (^) are thermal masses or soft SUSY masses, but are not due to the Higgs vev. 
Neglecting for the moment the mass effects, one hnds in the SM (with uh Higgs doublets) 

Qera « + Sur - Zdn - ~ ‘^riuH = 0 (23) 


In the MSSM, the numerical coefficients change slightly, assuming msusY ~ Tc^ because there 
are now two Higges and spartners as well as partners in the soup. The net effect is to change the 
number of Higges uh from one to two in (^). 

There are now seven homogeneous equations for the 10 unknown chemical potentials. If one hxes 
the asymetries in the three {H/3 — Lj}s, all the chemical potentials are determined. 

The baryon asymmetry carried by the quarks is 

^2col X)gen(2 X ('94') 

(X ^ > 


where I have used equations (O) and neglected 0{firn?) corrections. Quarks of different families 


have the same chemical potential so quark mass effects in equation (^) cannot be used to preserve 
an asymmetry. Equation ( ^^ can be written using equations (0), (ED and (^D 


ub = 


2A + Ann 
66 -I- 13n 


H 


^BI3-L, + 


18967r2 


{nie-Lr + 


(25) 


where uh is the number of Higgs doublets. I have included here the lepton mass effects—specihcally 
the r yukawa Hr = mT-/(175 GeV), because it is largest. These small corrections, which are present 
below 0 and above[E^ the EPT, will be discussed in section four. It is clear that if H — L = 0, 
one needs very large lepton flavour asymmetries (n^.-L./s ~ 10“^, where s is the entropy density) 


to preserve the baryon asymmetry via lepton mass differences. (As noted in |E^, slepton mass 
differences could be more effective.) 


7 









2.3 bounds on B — L violating coupling constants 


Equation (pSf ) says, as expected, that a baryon asymmetry will remain in the plasma if an asymmetry 
remains in any one of the {-B/3 — Lj}. Neglecting the lepton mass effects, it is clear that we need an 
asymmetry in at least one of the {-B/3 — Li} to preserve the asymmetry in the quarks. This means 
that interactions violating one of the {S/3 —Lj} need to be out of equilibrium. Suppose, for instance 
that it is S/3 — Le that is not washed out. Then one requires 


S/3-Ls 


< H 


(26) 


for Tc < T < Te^ ~ 100 TeV. 

I can translate (^) into a constraint on coupling constants using the rate estimates from section 

2 . 1 . 

For majorana neutrino masses that are present in the thermal soup above the EPT as a mass, 
requiring T from equations 


or (|I5|) to be less than S at T = 100 GeV gives 


m,_ 


< 10 - 100 keV 


(27) 


(The weaker bound corresponds to the scattering interaction rate). 

In the case where rriiy ~ above the EPT there is a lepton number violating scattering 

process mediated by the dimension 5 operator {iiHH)/M. As discussed in the introduction, /H 
increases with the temperature for D > 4 operators, so the best bound is at the highest possible 

0 - 


temperature. This is 




100 TeV 


T ~ 100 TeV gives the bound 


m. 


Requiring equation ([T^) to be out of equilibrium at 
, ^ keV . (28) 


These are bounds on a neutrino mass that violates electron lepton number. So in the basis where 
the charged lepton mass matrix is diagonal the majorana neutrino mass matrix entries m®* < 10 keV 
for z = e, /i, r. 

The i?-parity violating masses ejh 2 ^i are slightly trickier ||29|, because one can set them to 
zero by redehning the higgsino hi. The “problem” here is that the Hi superfield has the same gauge 
quantum numbers as the lepton superfields L*. In a lepton number conserving theory, Lj are distinct 
from Hi because they carry lepton number. However if the zth lepton number is not conserved, then 
there is no unique dehnition of which linear combination of Li and Hi is “the Higgs”. 

For setting bounds on lepton number violation, it should not matter which basis we calculate in. 
Whether an asymmetry in, for instance L^, survives in the plasma is a physical question that should 
not be basis dependent. It is possible to construct basis independent “invariants” that parametrise the 
amount of R parity violation between different coupling constants (analogous to Jarlskog invariants 
for CP violation), and to set bounds on the invariants |^, However, one can get the right 

answer with a less formal approach. 

If we rotate the superpotential term ejiJ 2 -hi into HH 1 H 2 by redefining Hi. 


H'i = 


y/JjFTB 


{fiHi + ejLj) 


(29) 


then we generate new trilinears: 




LiQpDl 


hjB, 


LiLjE 


C 

k 


(30) 


If the soft masses BhHiH 2 + BiH 2 Li are parrallel (in Hi, Li space) with (/i,ej), (he. Bi/Bu = 
ei//i), then they are simultaneously rotated away. If in this basis the soft masses mixing Hi and the 





















Li are also zero, then there is no lepton nnmber violation in the masses. In this case, there are two 
possible bonnds one can set on ep one conld require that the rate associated with the mass e, be 
out of equilibrium, or that the trilinear interaction < H. The hrst bound would 

imply e* < 10 keV, the second <10 ^/i ~ MeV. So it matters which is the right one to impose. It 
turns out that the second, weaker constraint is the correct one. One can see this in a one generation 
model; the stronger rate due to the mass term chooses the Higgs to be the direction of equation (|^) , 
so the lepton is the orthogonal direction, and lepton number is violated by the X'LQD^ term in this 
basis. One can check that in the basis the trilinear A can be written in terms of the coupling 
constants of the original basis as — /xA . Equation (l37|) gives the bound 

{eih^^ — fiX < 10“’^ for one i (31) 


or e* < MeV for /x ~ 100 GeV. 

It is unlikely in realistic models that (/x, e^) should be exactly parrallel to {BH,Bi), or that the 
soft masses should be zero in this basis. Small lepton number violating masses parametrising 

the misalignment should remain. For simplicity, let me neglect the and just consider e* and 

Bi. The rate for these masses can be estimated as a mass correction times a gauge interaction rate, 
like the estimate for the neutrino majorana mass, that is 


T ~ min 



/i2 + C/2T2’ 


{BH + 9^T^)^j 


X 10-V^ 


(32) 


where in this formula e, {Bi) is taken in the basis where Bi (e*) is zero. The denominator contains 
a crude attempt to include thermal mass effects. To see why equation (^^ is the right bound, 
imagine sitting in the early Universe above the EPT. As the temperature drops, more and more 
renormalisable interactions come into equilibrium. At some point the larger of the two mass terms 
{Bh)'^ + (i?i)2 and 9 ? + comes into equilibrium, and thereby chooses the Higgs direction— i.e. if 

Y^ {Bh)'^ + {Bi)‘^ > jj? + el the Higgs will be combination B^H + BiLi. then in this basis chosen by 
the interactions, one needs the lepton number violating rate associated with the mass e* to be out of 
equilibrium, which is equation (P^). 

The bound (p^ can be expressed in a more basis independent way as: 


Bnej - nBj 


< 2 X 10"^ 


(33) 


or, in the basis where Hcj is zero, ej//i < 2 x 10“^. So the cosmological bounds require that (/x, e*) 
be more aligned with {Bh, Bi) than with the yukawas. 

In the supersymmetric Standard Model, lepton flavours are unlikely to be separately conserved. 
The slepton mass matrix is unlikely to be exactly diagonal in the lepton mass eigenstate basis. If 
lepton flavour violating processes due to flavour off-diagonal slepton masses are in thermal equilib¬ 
rium, then all H — L violating processes need to be out of equilibrium, not just interactions violating 
on of the H/3 — L,, i.e. the bounds (|^) , (|37D and (p^) would apply to all H — L violating coupling 
constants, and not just those involving one lepton generation. So if one wants to use the “flavour 
loophole”, and only impose the bounds on one generation, the flavour changing soft masses must be 
out of equilibrium. 

The rate associated with a flavour off-diagonal soft mass can be estimated analogously to the 
previous discussion about ejif 2 -hi. There are three interactions that choose a basis in lepton flavour 
space: the lepton yukawa hl^, the A term and the soft masses mly If I want to separately 
conserve lepton flavours in the thermal soup, then all the interactions in equilibrium must agree on 
the basis in lepton flavour space. Suppose I start in the basis where is diagonal, because this is 
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a familiar basis from phenomenology, and I put a hat on the coupling constants in this basis. There 
will be off-diagonal terms [rh‘^]ij,A^J, i ^ j in this basis. Since the interaction rate associated with 
is larger at T ~ 100 GeV than the rate for h^, this is probably not a sensible basis. I should rotate to 
the basis where is diagonal, because the stronger interaction will choose the lepton flavour basis. 
In this new basis, there will be flavour off-diagonal yukawas ~ /m\ {i ^ j, and no sums on 

repeated indices). These yukawas need to be out of equilibrium: 


T ~ 10“^ 



2 

T <H 


(34) 


If one wants to preserve B/3 — L^, this gives the bound < 5 x 10“^. Preserving an different 

lepton flavour would give a stronger bound (the yukawa is larger) on different elements of rh^. 

Now consider the A term. Suppose I estimate the associated rate from equation (^) to be 


This is equivalent to assuming that the A term mediates a decay of one scalar into the other two. It is 
not clear that this is kinematically allowed, but I use this estimate anyway (I estimate the scattering 
rate r(00 —>■ 07) ~ 10 ~^g^A^/T, which gives a bound on A that is about an order of magnitude 
weaker.). In the basis where is diagonal, I need 


(z^j) (36) 

for whichever lepton flavour i one wishes to conserve. This gives A^^ < 10“® GeV. This is not a 
particularity strict bound; recall that I have absorbed the yukawa into A (see equation (||)). Similar 
bounds apply to the /^-violating A terms of equaltion (0). 

Finally I need to calculate a bound on B — L violating trilinear interactions. There are no basis 
confusions for these interactions, so to conserve, for instance, B/3 — L^, I need all the trilinears that 
violate Le and B to satisfy 

25T^ 

10“^A^T < FT ~- (37) 

rripi 

This means that < 10“^ for all p, g, r, Ag^^ < 10“^ for all p, g, and X^jk, Xjke < 10“^ for j, k ^ e. 


2.4 summary of bounds on R violating interactions 

Suppose that I dehne the Higgs such that the three Bi are zero. Then to conserve B/3 — L^, I need 

Ce < 10 - 100 keV, < (60 - 90 MeV)^ (38) 

for ~ ~ i? ~ 100 GeV. 

I need all the trilinears and A terms that violate B/3 — Le to respectively satisfy 

A < 10“^ , H < 10 - 100 keV (39) 

I also need Le flavour violation mediated by the R conserving soft masses to be out of equilibrium, 
which will be the case if 

9 

Tfl 

—f < 5 X 10-2 , 4“ < 10 - 100 keV (40) 
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3 if the BAU was not there... 


I wrote in the introduction that I would assume that the BAU was present in the thermal soup above 
the EPT. This is of course one of the loopholes in the constraints that I am discussing; if the BAU 
is not there, then there are no constraints on B — L violation. So in this section I will relax this 
assumption. 

First I would like to work out over what temperature range the BAU needs to be “not there”. 
The lower end of this range is straightforward to End. The B + L violating electroweak effects drop 
out of equilibrium at or shortly after the electroweak phase transition. Once they are gone, the BAU 
is no longer at risk. So the baryon asymmetry can safely be present in the plasma at T <C T^. 

It can also be present at high temperatures, before the electron yukawa comes into equilibrium. 
The interaction rate associated with a dimensionless couplings constant A scales linearly with the 
temperature: T ~ A^T. Since H ~ T'^/rUpi, this means that yukawa interactions are out of thermal 
equilibrium at sufficiently high temperatures, and come into equilibrium as the temperature drops. 
Since the electron yukawa is the smallest coupling constant in the SM, it is the last to come into 
equilibrium at Te^ ~ 100 TeV. The carries hypercharge, and the Universe ought to be hypercharge 
neutral, so if there is an asymmetry in the en, and no interaction in equilibrium that can transfer 
this asymmetry to other particle species, then asymmetries must remain in other particles to ensure 
that Y = 0. It is easy to check from the relevant equations of chemical equilibrium that this means 
a baryon asymmetry will remain in the plasma. At T ^ Te^, there will be other yukawa interactions 
out of equilibrium, and additional symmetries that can protect the baryon asymmetry, but the 
singlet cr is the last particle species to come into chemical equilibrium, so can protect the baryon 

0 


asymmetry for the longest 


Any B — L violating operator not involving the cr can be in 


equiliibrium above without being able to eat the BAU. This was not fully realised in the earlier 
papers discussing bounds on non-renormalisable B — L violating operators, where strong bounds were 
set on iA > 4 operators by requiring that they be out of equilibrium up to Treheau or the temperature 
when the BAU was generated, or the temperature when the B + L violation came into equilibrium. 

Note that this whole argument assumes that there is an asymmetry stored in the cr. This is 
maybe not so easy to arrange—for instance the decay of heavy singlet (“right-handed”) majorana 
neutrinos generates an asymmetry in the lepton doublets, but not the charged singlets. However, 
for setting generic bonds, one must allow for the possibility that some baryogenesis mechanism does 
generate an or asymmetry, which can then protect the baryon asymmetry down to Tg^. 

To approximately determine Tg^, one can assume that the cr comes into chemical equilibrium 
when the decay of the Higgs into a left handed and a right-handed lepton is of order H [O: 


T ~ IQ-^hlT ~ H 


(41) 


which gives Tg^ ~ 100 TeV. See for a careful determination of Tg^, that includes the various decay 
and scattering interactions mediated by hg, and the thermal masses of the participating particles. 

So the baryon asymmetry is really only “at risk” over three decades in temperature: Tg < T < 100 
TeV. If the asymmetry is present in the thermal soup during this period, it needs to be an asymmetry 
in B — L, and interactions violating at least one of the B/3 — need to be out of equilibrium during 
this period. This gives the bounds of section 2. 

Alternatively, the BAU can be not in the soup between Tg < T < 100 TeV. The most obvious way 
to do this is to generate the baryon asymmetry at the electroweak phase transition (see, e.g [|], |^ for 
reviews of electroweak baryogenesis). But this is not the only possibility; if we turn out to be living 
in a corner of SUSY parameter space that does not allow electroweak baryogenesis, there are models 
where the baryon asymmetry is generated in the late decay of some particle, after the EPT . 

Another possibility is to generate the baryon asymmetry early, and then “hide it” between T, 




and Tg. One way to do this |^, is to add additional particles that carry baryon number to the SM 
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or the SSM, decouple them from the rest of the SM particles at T > and have them decay after 
the EPT, returning the baryon asymmetry that they carry to the plasma of SM particles. Another 
possibilityy is to store some hypercharge in a coherent hypercharge magnetic held for Tc < T < 




If there is a sufficiently large asymmetry in the eji, it can act as a source for the hypercharge 
magnetic held before and this held may survive until the EPT. 

Another somewhat related loophole is to take the up-quark yukawa to be considerably smaller 


than is usually assumed |^. If hu ^ 10“^ (m^ ^ 20 keV !), then an asymmetry in the urs could play 
the same role as the cr asymmetry does, but all the way down to the EPT. If hu did not come into 
equilibrium until after the EPT, and if the urs carried a primordial asymmetry, then this asymmetry 
will be present after the EPT, irrespective of how much B — L violation is in equilibrium above the 
EPT (providing of course that the B — L violating operators do not involve ur). 


4 other loopholes 

In this section, I list the remaining loopholes. 


4.1 lepton mass effects 

In models where ur-l = 0, lepton mass effects can protect the baryon asymmetry from SM B + L 
violation, provided that there are lepton flavour asymmetries [^. A small baryon asymmetry will 
remain below m and above the EPT, proportional to lepton mass differences and flavour 
asymmetries. Since I am here discussing what happens above the EPT, the numerical factors in 
equation (^Sf) are for this case, where the lepton masses are “thermal”. In this subsection, I neglect 


most of the numerical factors. 

The number density for particles in thermal equilibrium is 


n = 


2 g r p^dp 

TT^ J X 


(42) 


where the ± is for fermions or bosons, and the number density h for anti-particles has the sign in 
front of the chemical potential p flipped. Expanding n — h in p/T and assuming that both 

are small, gives equation (^). As previously noted, to correctly determine the 0{pim?) term, one 
must calculate the free energy, and take the derivative with respect to p. However, I am not here 


interested in the exact coefficient; see ^ for numerical factors. 

The reason that these {pm?) contributions are interesting is that they can preserve a baryon 
asymmetry in the presence of the B + L violating interactions, when B — L = 0, and there are flavour 
asymmetries amoung leptons (also when B + L and B violation are in equilibrium, and there are 
flavour asymmetries). They do not protect the BAU if there are interactions in equilibrium that take 
all three lepton chemical potentials to zero. 

So suppose that the B + L violating interactions are in equilibrium (equation ^), and that 
Ur-l = 0 : 


12qL-3Y,. 



3H = 0 


(43) 


This could be the case if the BAU was generated in the out-of-equilibrium decay of heavy particles 
from Si B — L conserving GUT. There are two types of constraint imposed on the plasma: the sum of 
the chemical potentials participating in each interaction must be zero, and certain charge densities 
{Y,B — L) must be zero. If one includes the 0{pm?) terms in the charge densities, one hnds that 
there is a small BAU in equilibrium. This is due to the non-zero chemical potentials and not equal 
masses for the different lepton generations. 
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From the condition of charge neutrality, equation 
in terms Ql and 


H ~ 


, I can solve for the Higgs chemical potential 
3gL - E. m + OK/T^)] Uql -EiiiX OK/T2)] 
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(This assumes SM particle content, and I have dropped Hmj terms). Substituted into 
using (|I^ one hnds 

nB oc qlT^ ~ 


(44) 
, and 

(45) 


7 ^ 0, if there are lepton flavour asymmetries {nL._L^ ^ 0), so the remaining baryon 
asymmetry (^) is proportional to lepton mass differences, and flavour asymmetries, see equation 
(|25|). It is clear from (p5D, where all the numerical factors are included, that one needs very large 
lepton flavour asymmetries (~ 10“^ — 10“^) to preserve a sufficiently large baryon asymmetry. An 
interesting possibility suggested in is to use the slepton mass differences in a SUSY model, 
because these could be much larger. 


4.2 modified cosmologies 

In this subsection I briefly review some modihcations of cosmology that allow one to escape the 
bounds discussed in section 2. 

One possibility is low T^eheat models; if the Universe after inflation reheats to a temperature 
Treheat < T^, then oue never has to worry about having B + L and B — L violating operators simulta¬ 
neously present. However, one does have to worry about how to generate the baryon asymmetry in 
such models. This is not neccesarily a problem; one can, for instance, generate an asymmetry from 
the decay of the inflaton, or in some cases via the Affieck-Dine mechanism [^. In fact, the reheat 
temperature after Affieck-Dine baryogensis is often low. 

Another possibility is to modify the expansion rate of the Universe. If the Universe expanded 
somewhat faster, the cr could protect the BAU all the way down to Tc. This can be accomplished 
by, for instance, adding a component to the energy density of the Universe that scales as 1/a « H 


(This is not as far-fetched as it might sound; a scalar held in the right shaped potential can behave 
this way.) 

It is also possible to “get rid of” the B + L violating electroweak interactions. For instance, if 
there is a SU(2) doublet with a large vev above the EFT, it would give a mass to the W, so the 
rate for B + L violation would be exponentially suppressed. This could occur in an an Affieck-Dine 
baryogenesis scenario, when the condensate carries a large baryon asymmetry 
does not decay until T ~ 100 GeV. 


39 , and therefore 


5 summary 

The argument leading to cosmological bounds on B — L violating interactions goes as follows. Assume 
that the baryon asymmetry was created before the electroweak phase transition as an asymmetry in 
at least one of the three H/3 — Lj, and that it is present in the thermal soup above the electroweak 
phase transition. The electroweak B + L violating interactions are operating in this plasma. If there 
are also interactions in equilibrium that violate all of the {B/3 — Li}, the baryon asymmetry will 
be washed out. Since it is supposed to be around today, this is undesirable. Therefore one requires 
that interactions violating one of the B/3 — Li be out of thermal equilibrium, so they cannot eat the 
baryon asymmetry. 
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This generically gives bounds of order A < 10“^ on “yukawa-type” coupling constants, and m ^ 
keV-100 MeV for different masses. The mass bounds can be slightly tricky to calculate. These 
constraints apply to interactions that violate B and one particular lepton flavour. 

There are various loopholes to this argument. One can avoid them completely by not having 
the BAU present in the plasma for 100 TeV > T > 100 GeV. This is approximately when it is 
at risk from the B + L and B — L violating interactions. The most obvious way to not have the 
baryon asymmetry present is to generate it at the electroweak phase transition. There are other 
models where the BAU is generated after the phase transition, or where it is generated at very high 
temperatures and hidden for 100 TeV > T > 100 GeV. 

It is also possible to modify the cosmological model to avoid these bounds. If the reheat tem¬ 
perature after inflation is below that of the electroweak phase transition, there is no opportunity to 
wash out the BAU. (One has fewer ways of making it though). If there electroweak B + L violating 
interactions are turned off above the phase transition, then an asymmetry in i? -|- L can survive. If 
the expansion rate of the Universe is increased, the slowest Standard Model interaction, which is 
the electron yukawa, might not come into equilibrium until after the phase transition. If there is 
an asymmetry in the e^i, this means that an asymmetry must remain in the baryons, because the 
Universe should be hypercharge neutral. 

It is possible to avoid the strong bounds on B violating operators, or the conclusion that the 
asymmetry must carry B — L, via lepton mass effects. If there are large lepton flavour asymmetries, 
and lepton flavours are conserved, then a small baryon asymmetry proportional to m^/T^ x nL^_L. 
will remain, even if ub-l = 0 or if there are B violating interactions in equilibrium. 
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